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Abstract 

Data are presented on the reaction pp — )• ujujtt^ at rest from the Crystal Barrel detector. 
These data identify a strong signal due to /2(1565) — )• ojoj. The relative production from 
initial pp states ^P2, ^-Pi and is well determined from ujuj decay angular correlations; 
P-state annihilation dominates strongly. A combined fit is made with data on pp — )• 3tt^ 
at rest, where /2(1565) — )• vr'^vr'' is observed. A Flatte formula is fitted to the /2(1565), 
including the s-dependence of decay widths to ujuj and pp. The data then determine the 
K-matrix mass, M = 1598 it ll{stat) it 9{syst) MeV. The decay width to 2tt is very small, 
of order 2% of the total width. 

We investigate the /2(1565) through its decays to uu in the process pp — )■ uutt^ at rest. Data 
on pp Stt" at rest are also used to estimate relative widths to vrvr and cou. 

The first clear identification of the /2(1565) was by May et al. in the Asterix experiment on 
pp — > 7r^7r+7r° at rest [1,2]. Subsequently, it has been identified by the Obelix collaboration in 
np 7i~7i~^TT~^ [3], and in Crystal Barrel data on pp — )■ 37r° at rest [4,5]. However, it has become 
clear that the /2(1565) also couples strongly to uu (and therefore probably to pp). Abele et al. 
[6] found a strong cusp in the Tin D-wave at the uu threshold in the 37r° data. This led to the 
identification of /2(1565) as the same 2"*" resonance as observed by GAMS [7] and VES [8] in uu 
just above threshold. Their results are listed under /2(1640) by the Particle Data Group (PDG) 
[9], who also list other sightings there and under /2(1565). The properties of this resonance have 
remained elusive. Here, we show that it is important to fit with a Flatte formula which includes 
the s-dependence of decays to vrvr, pp and uju; there is a sharp cusp at 1564 MeV due to the 
opening of the uu channel. 

In outline, the elements in the analysis are as follows. Firstly, uu correlations in the Crystal 
Barrel data on the uun^ channel determine well the relative amounts of annihilation from initial 
states ^5*0, ^Pi and ^P2- Secondly, earlier Crystal Barrel data on the 37r'' channel determine well 
the lower edge of the resonance. Both of these sets of data are readily fitted with a Breit-Wigner 
amplitude which includes the s-dependence of widths for decays to pp, ujuj and tttt. However, 
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the upper side of the resonance is obscured in Crystal Barrel data by the centrifugal barriers 
associated with the production reaction; there is an L = 2 barrier for production from ^5*0 and 
L = 1 barriers for production from '^Pi and '^P2. These barriers cut off the resonance on its 
upper side in pp data at rest. The data from GAMS and VES for the uu channel provide a 
good determination of the upper side of the resonance. We include the VES data of Ref. [9] 
into the analysis in order to provide this constraint. 

Firstly, we present technical aspects of data selection. The new uutt^ data reported here 
concern events where both u decay to 7r"'"7r~7r°, so the final state is 27r+27r~37r°. They were 
taken with the Crystal Barrel detector [10] using a p beam of 200 MeV/c stopping in liquid 
hydrogen. In early runs, the target was surrounded by two cylindrical multiwire chambers for 
triggering, followed by a cylindrical jet drift chamber measuring charged particles. In later runs, 
the multiwire chambers were replaced by a silicon vertex detector. The outermost layers of the 
jet drift chamber covered ~ 70% of the solid angle. Surrounding this drift chamber was a barrel 
of 1380 Csl crystals detecting photons with good resolution and high efficiency over a solid angle 
98% of An. A solenoidal magnet provided a field of 1.5T. 

The present data sample consists of 9.4M triggered 4-prong events. The trigger requires 4 
hits in the multiwire chambers (or Si vertex detector); events with 4 long tracks are selected by 
demanding a hit in the first 3 and last 4 layers of the jet drift chamber. Off-line analysis selects 
~ 2.4M events with four well reconstructed tracks and balancing charges. Gamma rays are 
selected by demanding a shower with energy larger than 10 MeV. Any energy deposits matching 
up with charged tracks are rejected. To eliminate the risk of shower energy being lost down 
the beam-pipe, events are rejected if the centroid of any shower lies in crystals immediately 
adjoining the entrance or exit beam-pipe. These cuts reduce the sample to ~ 1.3M. 

Events with 4 charged tracks and > 67 are then submitted to a 7C fit to 2n^27i~3n^; 11,679 
events survive with a confidence level above 0.1%. Almost all come from events containing only 
6 or 7 showers. In the latter case, one shower is interpreted as arising from a so-called 'spht-off ' 
in the Csl crystals. This is caused by nuclear interaction of a charged particle, which generates 
neutrons that undergo secondary interactions in nearby crystals; those low energy deposits close 
to charged tracks are discarded, but some neutrons convert far from parent tracks. 

Turning now to physics results, we consider first the branching ratio to ujoutt'^. The Monte 
Carlo simulation is unable to provide a precise efficiency for the reconstruction of charged tracks, 
due to the 'split-offs'. Accordingly, we normalise the cucutt" branching ratio to that for urjii^ . 
The latter is known from an earlier study of 77 events [11] and is (6.8 ± 0.5) x 10^^. Fig. 1(a) 
illustrates the mass spectrum of 7r"'"7r~7r° combinations. We estimate 930 ± 50 rjujn^ events and 
2000 ± 50 uu-K^. Correcting for the branching ratios BR[u 7r+7r°7r°] = 0.888 ± 0.007 and 
BR[r] 7r+7r°7r°] = 0.274 ± 0.026, we find a branching fraction 

BR^p ujujTT^] = (4.5 ± 0.7) X 10"^ (1) 

We next consider combinatorics and backgrounds. The kinematic fit to cocott^ has 12 different 
combinations of pions (and 180 permutations including individual photons). It is possible that 
in some cases the confidence level of an incorrect combination will be higher than that of the 
correct combination. To reduce this problem, two further cuts have been applied, based on 
Monte Carlo simulations. The first is that there should be two or less combinations above 5% 
confidence level. The second is that the ratio of the confidence level of the best combination 
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Figure 1: (a) The 7r"'"7r~7r° mass spectrum from events fitting the 27r"'"27r~37r'' hypothesis with 
confidence level above 0.1%; (b) the mass distribution of the second best cu combination. 
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to that of the second best should be above 0.6. With these cuts, the Monte Carlo simulation 
estimates that wrong combinations are less than 16%. They occur mostly amongst similar 
geometries, which the simulations show to have small effects on the fitted physics. 

A possible source of background to uun^ is the ujn~^n~2n^ channel. To investigate this, we 
form all possible combinations of 7r'^7r~7r°, select the best ou and then plot the invariant mass 
distribution of the second best combination. This is shown in Fig. 1(b). There is a dip at the co 
mass, due to selection of the best one. The second best combination peaks at the u, but there 
is a background which, extrapolated under the u, amounts to ~ 4%. A cut is applied such that 
the invariant mass of the second best u combination lies between 752 and 812 MeV, so as to 
reject obvious background events. Finally a kinematic fit is made to ujoutt^, and events fitting 
with confidence level above 5% are accepted. There are eventually 1346 accepted events and 
4677 from the Monte Carlo simulation. 

Fig. 2 (a) shows the Dalitz plot for data and Figs. 2(c) and (d) the mass projections for uir 
and combinations. Full histograms show the result of the maximum likelihood fit described 
below; dashed histograms show distributions derived from phase space after acceptance cuts. 
The Dalitz plot from the fit is shown in Fig. 2(b). There is a clear peaking of events towards 
the top right-hand edge of the Dalitz plot. This could be due to /2(1565) — ?> uu, /o(1500) — )■ couj 
or 6i(1235)a;. Distinctions may be made using angular correlations for decay of the two u. 

In addition, data are also available for the ojujti^ final state from 87 data, where both u — )■ n^j. 
The Dalitz plot for 740 events from annihilation in hydrogen gas is shown in Fig. 2(e). The 
/2(1565) appears prominently near threshold; it is enhanced in P-state annihilation in gas. These 
data have been shown to give amplitudes consistent with the analysis presented here. However, 
we shall not present their analysis in full. The reason is that much of the spin information about 
cjuj correlations is lost in the polarisations carried by decay photons from the a; decays. In this 
respect, charged decays of the uj are greatly superior, since no spin information is lost. 

The primary processes which we consider in the amplitude analysis are: 
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Here a denotes a slowly varying 0+ component decaying to lulu, parametrised as a constant 
amplitude. We have also tried contributions for production of lulu with — 0~, 1~ or 2~ 
resonances having widths set to 250 MeV, but find them to be consistent with zero within 
errors. 

For spin resonances, reactions (3)-(5), the amplitude is proportional to ua-Ub, where n = 
Pi Ap2 and pi^2 are 3-momenta of 7r+ and 7r~ from decays of each lua^b in their rest frames. For 
spin 2, amplitudes are given by complicated tensor expressions [12] in terms of ua, ub and the 
decay angles 9i of each resonance to lulu in its rest frame. Fig. 3 illustrates the projection of 
decay probabilities on to cos0, where (J) is the angle between the vectors ua and ub normal to 
the u decay planes. For spin 0, the distribution is simply cos^ 0; for spin 2, it is much fiatter in 
terms of 0, but in addition has a distinctive correlation between 0, 9i and 62 which is included 
in the amplitude analysis. The data clearly demand a mixture of 0"*" and 2+. 
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Figure 2: Dalitz plots for (a) data, (b) fit. Projections on to (c) M^(aj7r) and M^(a;a;); points 
with errors are data, full histograms show the fit and dashed histograms show phase space after 
acceptance cuts; (e) the Dahtz plot for 87 events pp — >■ uun^ in hydrogen gas. 
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As a check, we have examined four slices of uu mass and determined the intensities of the 
0"*" signal and 2+ originating from ^P2-, ^Pi and ^5*0 in each slice. The 2"*" signal peaks strongly 
near the ujuj threshold, while the O"*" signal extends fairly uniformly across the whole ujuj mass 
range. It is therefore clear that the /2(1565) accounts for a large proportion of the observed 
events. This is also evident from Fig. 2(e). 

Some further simplifications in the analysis are possible. We find that there is little distinction 
between cttt and /o(1770)7r final states, because of the limited ujuj mass range. We therefore drop 
/o(1770). Secondly, there is little distinction between the /o(1500)7r channel and 6i(1235)c<;. The 
reason is that both processes involve mostly orbital angular momentum L = in the production 
reaction and both enhance the upper right-hand edge of the Dalitz plot. [In fitting 6i(1235), 
we allow for decays to utt with both i = and 2 and take the small (0.29) D/S ratio of 
amplitudes from the PDG value]. The 6i (1235)0; channel gives a somewhat better log likelihood 
than /o(1500). However, if both processes are introduced, they fit to large magnitudes with an 
unreasonably large destructive interference between them. It is therefore necessary to make a 
choice between them. We now argue from other data that the contribution from /o(1500) is 
small. 

Firstly, in amphtude analyses of pp 37r° at rest [4-6], there is no indication of any cusp in 
the /o(1500) amphtude at the uu threshold. Secondly, data on J/'^ 7(47r) require /o(1500) 
decaying dominantly into aa rather than pp [13]; at the quark level, coupling to pp is three 
times that to ujuj, making the ujuj decay weak. Thirdly, analysis of Crystal Barrel data at rest 
to Stt shows a large signal for /o(1500) — )■ aa, but little or none for decay to pp [14]. From these 
sources, we estimate that the branching ratio of /o(1500) to uju is less than 4% of that to tttt. 
The branching fraction for pp — )■ /o(1500)7r followed by the decay /o(1500) — >■ tttt is known from 
Sn^ data to be ~ 2.45 x 10~^ at rest [6]; hence the branching fraction for its decay via uu is 
below 10~^. Using equn. (1), /o(1500) — >■ uu must be less than 2.1% in intensity in present data 
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and we drop it from the analysis. Including it with an intensity 3 times this limit produces only 
a small interference with the 6i(1235)a; final state and has little effect on conclusions concerning 
/2(1565). 

Contributions to the fit are shown in Table 1. Those from cttt and bi{1235)u} are lumped 
together because of interferences which are irrelevant to the present discussion. It is obvious 
that P-state annihilation dominates production of the /2(1565). The ratio of intensities from 
and ^Pi is close to their multiplicity ratio 5/3. Errors are mostly due to statistics, but also cover 
variations of centrifugal barriers and form factors discussed below. The first three components 
of Table 1 are very stable in intensity, whatever changes are made to other components; they 
are also insensitive to the choice between a Breit-Wigner amphtude or a Flatte form in fitting 
/2(1565) and to its precise mass and width. This insensitivity refiects the fact that the ojoj 
correlations are able to identify the spin dependence distinctively. 

Channel contribution (%) 
^-^0 /2(1565)7r (4.2 ±0.2) 

^^1 ^ /2(1565)7r (19.5 ±2.8) 

^^2 ^ /2(1565)7r (30.9 ±3.3) 

^ 6i(1235)a; + a7r (45.4 ± 4.2) 



Table 1: Percentage contributions to a;a;7r° data. 

We discuss now the formula used to fit production of the /2(1565): 

^ ABl{p)V{p) exp(-ag^) 
/1565 M-' - s - m{s) - iMFtot' ^ ' 

Ttot = ^27r + 9lPl + 92P2- (8) 

In the numerator of equn. (7), V{p) is the Vandermeulen form factor for the production process 
pp — >■ /2(1565)7r [15]; p is the centre of mass momentum of /2(1565). The factor Bl{p) is the 
standard Blatt-Weisskopf centrifugal barrier for production with orbital angular momentum L. 
Elsewhere, we find a radius for the centrifugal barrier R = 0.5 to 1.0 fm; here, we extend the 
possible range to 1.3 fm, to accomodate the possibility that /2(1565) is an cucu 'molecule' of large 
radius. Fits to data are insensitive to R, but it will play a role in unfolding the fine-shape of the 
resonance. The exponential is a form factor with a = 1.5 GeV~^ for decay of the resonance to 
uu, each of momentum q. The factor A is a complex coupling constant for ^5*0 production; for 
P-state production it may be taken to be real, since only /2(1565) contributes. Other processes 
are described by equations analogous to (7), except that -Bl(p) = 1 and T-tot is constant. 

In equn. (8), the 27r width is taken to be constant. We shall find that it is very small. Care 
is therefore necessary over the s-dependence of the widths for decays to pp and ujuj. The factors 
Pi^2 of equn. (8) arc phase space factors for pp (channel 1) and uuj (channel 2). For the ujuj 
channel, p2 = 2gexp(— 2Q;g^)/y^. For the pp channel, the pi{s) is evaluated numerically from 
equn. (40) of Bugg, Sarantsev and Zou [16]. It is illustrated in Fig. 4. It has been parametrised 
by the following empirical expression: 

-3.909 + 10. 57Ls - 1.81s2 
^^^^^ ~ 48.47[1 ± exp(ll. 353(1.063 - s) + 5^(4.572 - 0.826s))] ' ^ ' 
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with s is GeV. The quark model predicts gi — 3g2, corresponding to the three charge states of 
pp and one for cucu. We explore other values of the ratio gi/g2 — 1 to 4, because of the possibility 

that /2(1565) is a 'molecule'. 

In equn. (7), m(s) is an important dispersive correction to the mass evaluated from the 
subtracted dispersion relation: 

mis) = - / ""^'"'^''^ (10) 

Here, Ttot is the total width appearing in cqun. (8). This dispersive correction makes the 
amplitude fully analytic. It has been checked by evaluating the standard dispersion relation for 
the real part of the amplitude in terms of an integral over its imaginary part; this relation is 
accurately satisfied over the whole range of s relevant here. 

Our objective is now to make a combined fit to the uun^ data and 37r° to determine M, r2-K, 
gi and g2- We also use VES ucu data to constrain the width of the uu peak at half-height. In 
practice, what we find is as follows. Firstly, the uojit^ data depend fairly weakly on the mass 
of the resonance and its full width. Their main role is to determine the magnitudes of the ^5*0, 
^Pi and ^P2 components. Secondly, the 37r° and uun^ data cannot be fitted simultaneously 
with a simple Breit-Wigner amplitude of constant width; this would require differences between 
channels for mass and width outside errors. They are, however, readily fitted together by equns. 
(7) and (8). 

The fit to VES data and 37r° determines well the mass and the full width of the resonance 
Ttot. The 37r° data define the lower side of the resonance and VES data the upper side. The 
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fit to 37r° uses the same ingredients as Abele et al. [6]. It is sensitive to the ^Sq component 
of production of /2(1565), since this component interferes with other strong amphtudes. It is 
insensitive to the magnitudes of the and "'Pi amphtudes for producing /2(1565), since they 
interfere only with the weak /2(1270)7r channel; the intensities of '^P2, ^Pi — )■ /2(1565) are small 
(~ 1.2 and 0.7% respectively). [They are smaller with respect to ^Sq than in Table 1 because 
the TTTT data extend to lower masses than ujuj and the centrifugal barriers therefore create less 
suppression.] 

The value of M optimises at 1598 MeV, with a statistical error of ±11 MeV. Systematic 
variations with R and gi/g2 are ±9 MeV. The dispersive correction m{s) plays a vital role in 
this stability. As M is varied, the subtraction point in the dispersion integral alters and m{s) 
moves up and down bodily at all s; the data will not tolerate much movement of this term, 
which is illustrated in Fig. 4(b). Without m{s), the quality of the combined fit deteriorates 
sharply. A detail is that the narrow spike in m(s) between 1.56 and 1.6 GeV originates from 
the rapid opening of the uu channel. 

The full width at half maximum of the summed cross section for tttt, ujuj and pp channels is 
determined reliably as 220 ± 15 MeV. For gi/g2 = 3, this corresponds to gi — 435 ± 30 MeV, 
with a systematic error from variations of R and gi/g2 which is negligible compared with the 
statistical error. However, the split between Tpp and F^^ is uncertain. 

The value of F27r is very small: 2.4 MeV. This small value follows from the very small intensity 
for production of /2(1565) in pp — >■ 37r° from ^5*0. This is (0.50±0.05)% of the 37r° channel; using 
the branching fraction (6.2 ± 1.0) x 10~^ for pp — )■ 37r° [17] and allowing for 3 charge states, the 
resulting branching ratio for production of /2(1565) in vrvr decays is (0.93±0.18) x 10"^. The rate 
for production of /2(1565) from ^5*0 and decay to uu is somewhat higher: 0.047 x 4.5 x 10~^ = 
2.2 X lO^''. However, this latter figure hides the fact that production of /2(1565) from ^5*0 
is highly suppressed by the centrifugal barrier at the high masses available for ujuj decay. To 
quantify this, we unfold the effect of the centrifugal barrier, and also the obvious suppression 
due to the limited phase space available in pp annihilation; the latter suppression factor is p/M, 
where p is the centre of mass momentum of the /2(1565) and M its mass. For gi/g2 — 3, we 
find for F27r, F^^^^ and Tpp, averaged over the line-shape of the resonance: 

f 2. : r^o, : r,, = 0.06i°:^2 : 1.23 ± 0.6 : 3. (11) 

The errors on r27r express a factor 3 uncertainty in r27r/rpp arising from uncertainties in the 
centrifugal barrier. There is also a factor 1.5 uncertainty in F^ju} from this source. Despite these 
sizeable errors, what is clear is that the 27r width is very small. This explains why the /2(1565) 
has not been observed in tttt — )■ tttt, e.g. in the CERN-Munich experiment [18]. 

In equn. (11), the ratio of pp widths to ujuj is of course dependent on gi/g2- What is accurately 
determined is the full width. The value of gi/g2 cannot be significantly less than 2, otherwise 
the TTTT channel is too strong at low masses and the fit to the Sn^ data deteriorates sharply. A 
small detail is that, in equn. (11), Tojoj/^pp is just above the value 3 assumed for gi/g2- The 
reason is that the finite width of the two p suppresses the pp signal slightly in the mass range 
of the /2(1565). 

To clarify these results, we display on Fig. 5(a) the shape of the resonance in 2n, ujuj and pp 
channels of production from ^5*0. The 27r result is what is fitted to 37r° data (scaled by a factor 
5 for purposes of display). It is strongly asymmetric, faUing to half- height at 1.59 GeV because 
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Figure 5: The intensity distribution for /2(1565) — > tttt (full curve), uu (dashed) and pp (dotted) 
for production in pp annihilation from ^Sq; (b) as (a) for the /2(1565) uninhibited by the 
centrifugal barrier effects or by phase space for production in pp annihilation. 



of the rapidly opening ujuj and pp channels. The ojoj result is likewise what is fitted to ojojtt'^ 
data for ^5*0 annihilation. It peaks at 1.59 GeV, but is rapidly attenuated by the centrifugal 
barrier. The pp channel peaks at the uu threshold; it falls on the lower side because of the 
falling phase space and on the upper side because of the centrifugal barrier and the opening of 
the uu channel. 

Fig. 5(b) shows our estimate of corresponding results for an /2(1565) in isolation, uninhibited 
by the phase space for production in pp annihilation or by the centrifugal barrier there. The 
results are subject to a factor 3 uncertainty of scale and some uncertainty of shape in unfolding 
the effect of the centrifugal barrier. Nonetheless, they illustrate one unavoidable feature. There 
is a cusp in the tttt channel at the (jjcu threshold, clearly visible in Fig. 5(b). The pp and cucu 
channels peak at ~ 1.63 and 1.66 GeV respectively. A Breit-Wigner amplitude of constant 
width is a poor approximation. 

In summary, we find that /2(1565) is produced in the pp — )■ uun^ channel with well determined 
branching fractions from ^5*0, ^Pi and ^P2; P-statc annihilation dominates strongly. It is well 
fitted by a resonance with strongly s-dependent widths in pp and cucu channels and a mass 
of 1598 ± ll{stat) ± 9{syst) MeV. Unfolding the effects of the uncertain centrifugal barrier is 
problematical, but the 27r width is certainly very small, of order 2% of the total. 

The /2(1565) makes a natural candidate for the radial excitation of /2(1270). However, the 
fact that its binding energy in pp and uu channels is small implies there must be some 'molecular' 
tail to its wave function, analogous to the long range tail of the deuteron. 
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